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Abstract 

Secondary instabilities of Faraday waves show three regimes: (1) As seen pre- 
viously, low-viscosity {v) fluids destabilize first into squares. At higher driv- 
ing accelerations a, squares show low-frequency modulations corresponding 
to the motion of phase defects, while theory predicts a stationary transverse 
amplitude modulation (TAM). (2) High-z^ fluids destabilize flrst to stripes. 
Stripes then show an oscillatory TAM whose frequency is incommensurate 
with the driving frequency. At higher o, the TAM undergoes a phase insta- 
bility. At still higher a, edge dislocations form and fluid droplets are ejected. 
(3) Intermediate-z^ fluids show a complex coexistence of squares and stripes, 
as well as stationary and oscillatory TAM instabilities of the stripes. 

47.20.-k, 47.35.+i 
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I. INTRODUCTION 



Over the past fifteen years, tlie problem of pattern formation and dynamics in non- 
equilibrium, spatially extended systems has attracted much attention. Although the most 
thorough and precise work has been on Rayleigh-Benard Convection (RBC) and Taylor- 
Vortex flow, another hydrodynamic system, the Faraday instability [0] , shows great promise. 
Here, surface waves are excited on a fluid by vertically vibrating its container. The instability 
is due to parametric resonance: in the reference frame moving with the container, gravity 
is modulated periodically at a frequency uo = 27r/, leading to surface waves vibrating at 
ujr = uj/2. Among the allures of the Faraday instability are its short time scales, large 
aspect ratios, and small number of material parameters. 

Large aspect ratios are important, for distant boundaries can still strongly influence 
pattern formation p[. A figure of merit is the aspect ratio T = L/^, where L is the system size 
and ^ the spatial correlation length. In many experiments, ~ A, the pattern wavelength. 
Aspect ratios of order 100 have been achieved recently in RBC and in chemical reactions 
A number of experiments on the Faraday instability using low- viscosity fluids have 
been conducted in containers with geometrical aspect ratios F' = L/A ~ 100 



However, as pointed out by Edwards and Fauve |jTO|, the damping length (f» ^) of waves 
in low-viscosity fluids such as water and butanol can be much larger than A (one can have 
F' ^ 100 while F ^ 1), implying that patterns in small-F experiments should be sensitive 
to details of the boundary conditions well above threshold, even when F' is large. Indeed, 



they are |11 



In this letter, we report observations of secondary instabilities for varying viscosities. We 
focus on high-viscosity fluids where ^ ~ A, in a box of moderate aspect ratio (F ranges up 
to 12, depending on uj). At high viscosities, the primary instability gives a stripe pattern, 
in contrast with the squares observed for low viscosities [0,|1^]. The sequence of secondary 
instabilities of stripe patterns is quite different from that observed for squares. At moder- 
ate viscosities, the competition between square and stripe states makes the situation more 
complex still. 

In low-viscosity fluids, squares destabilize via low-frequency modulations that Ezersky 
et al. interpreted as transverse amplitude modulations (TAM) of one component of the two 
sets of standing waves making up a square pattern [^]. Amplitude equations for standing 
waves do indeed predict a secondary TAM instability |]5| JT3|Jl^ . A similar analysis was 



assumed to hold for square patterns, but Douady [|15l showed by stroboscopic methods 
that the supposed TAM instability was actually a set of kink defects separating two sets 
of waves 180° out of phase, an observation that we have confirmed. Here, we show that 
when the primary instability is a stripe state, the secondary instability is an oscillatory 
TAM instability. A stationary TAM instability of very small amplitude is, however, seen in 
moderate-viscosity fluids. 

The above concerns control-parameter ramps, where the acceleration of the fluid con- 
tainer is slowly increased from zero. Quite different phenomena, such as converging target 
patterns and spirals, are observed when the control parameter suddenly jumps from zero to 
a large value |0|,|l6|,|l3 . 
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II. DESCRIPTION OF THE EXPERIMENTAL APPARATUS 



Our apparatus has been described at length in a previous study of acceleration thresholds 
and onset wavelengths [0. That study showed that under the proper conditions, we could 



predict quantitatively thresholds and onset wavelengths using a theory that assumed infinite 



lateral extent, infinite depth, but finite viscosity [T^. 

We used a square sample container of side L = 12.0 cm and depth 1.0 cm. The fiuid 
temperature could be controlled between and 60 °C ^0[, with a stability of ±0.001°C. 
Temperature variation due to vertical gradients was 0.5 °C in the stationary fiuid but was 
less than 0.05 °C once surface waves appeared. Our control parameters were the driving 
frequency u and the acceleration a, which was measured directly by an attached accelerom- 
eter. The acceleration could be varied from to 17 g's and was regulated to 1%. A spatial 
homogeneity of < 1% variation was achieved by careful leveling and dynamic balancing. 
"Soft" boundary conditions were created by using sloping sidewalls, which greatly reduced 
wave refiections from the boundaries. 

We used three different working fluids: a vacuum-pump oil made of paraffins (Fisher 
Scientific, Brand 19), propylene glycol, and glycerine. The properties of the pump oil were 
stable over days and allowed repeatable threshold measurements. Propylene glycol and 
glycerine could be doped with paint pigment, as described below, and allowed us to scan 
a wide range of viscosities (0.2-5 cm^/sec). Polar fluids, however, have much less stable 
material properties: their surface tension decreases upon contamination by small amounts 
of impurities, and their viscosity decreases if small amounts of water are absorbed. The 
sensitivity of the surface tension to impurities implied that different runs for glycerine and 
propylene glycol samples gave values of the surface tension that differed by as much as 
30%. In data quoted below, the material parameters were determined by flts to acceleration 
thresholds and onset wavelengths. These values were stable to about 5% over the course of 
a day. Again, the oil samples did not display such variability. All phenomena observed in 
the polar fluids were also present in the oil. 

We used three methods for visualizing the fluid's surface. Because the waves are very 
steep (surface inclinations can exceed 90°; see Fig. |I|f), methods based upon specular reflec- 
tion, such as the shadowgraph technique, cannot give much information about the detailed 
structure of the waves. (1) We doped the fluid with paint-pigment particles in order to turn 
the fluid surface into a diffuse light source [^. We added 2% (by weight) polymer-stabilized 
titanium dioxide particles to the glycerine and propylene glycol (the particles could not be 
suspended in the pump oil). We then illuminated the fluid surface using a stroboscope 
whose light was directed from the side of the container at an angle that minimized specular 
reflection into the CCD camera. (See Fig. |I|.) (2) We projected a known shape (straight 
edge) onto the pigment-doped fluid surface. From its apparent form, we could infer the 
surface proflle. (3) We suspended a ring of LEDs above the container []TB[. The LEDs were 
strobed with millisecond pulses and gave the positions of wave peaks and troughs. This 
method showed that all of the phenomena observed in glycerine and propylene glycol (with 
or without the paint particles) were also present in the oil. 
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III. EXPERIMENTS 



Figs. |I| and |^ summarize the patterns seen in high- viscosity fluids when the acceleration is 
increased slowly from zero. Below a critical acceleration Od, the surface remained flat. Just 
above a^i, for e = {a — ad) /ad ^ 10~^, we observed complicated, slowly shifting patterns 
[Fig. |l|a). Because the bifurcation away from a flat interface is (nearly) supercritical p2 



the correlation length is large when e ^ 0. (One expects ^ ~ ^okl~^^^5 where is the 
correlation length far from the transition.) In Fig. [l|a, ^ ~ (F 1.) The effects of spatial 
gradients in e are important in this region, further complicating the observed patterns. 

At higher a, the surface waves spontaneously order into standing waves of stripes when 
F ^ 3. (The F- value where stripes order depends on the boundary conditions.) (Fig. |I|b). 
In different runs, the stripes were parallel to one or the other of the sides of the square 
container. 

Above a second acceleration threshold, ac2 (the value of e depends greatly on uj), we 
observed a supercritical bifurcation to an oscillatory TAM instability of the stripe state. 
(Fig. |l|c). In one side of the cell, left-going traveling TAM waves dominate; in the other 
side, right-going TAM waves dominate. In the center, the modulations formed a standing 
wave. The frequency ujtam of the standing wave was incommensurate with the driving 
frequency /. (See Fig. |^.) The TAM frequency depends approximately linearly on /, but 
systematic deviations are clearly seen. The solid line in Fig. ^ shows a quadratic fit, which 
needs theoretical justification. The wavelength of the TAM was (3.30 ±0.35) A over the range 
of frequencies studied; the ratio uj/ojtam varied from 5.15 to 5.39 {i.e., ojr/ojTAM ~ 2.58 to 
2.70). 

At a third acceleration threshold, Ocs, we observed the spontaneous formation of de- 
fects. The line in Fig. |^ represents the acceleration where isolated single defects — edge 
dislocations — first appeared spontaneously. (Fig. |I|e.) They moved to the cell edge by a 
combination of climb and glide. Above Ocs, the number of defects increased rapidly, resulting 
in very complicated dynamics. (Fig. |l|f.) We have not analyzed this regime, but it might 
be useful to connect it with studies of defect-mediated turbulence 0. In the defect regime, 
droplets of fluid were spontaneously ejected upwards. Droplet emission was correlated with 
pattern defects, and their frequency rapidly increased with e. 

Between ac2 and acs, there is clear evidence for another threshold where the TAM mod- 
ulations dephase (Fig. |l|d.) However, the dynamics of this regime are complicated, and our 
cell was too small to study this tertiary instability quantitatively. 

We also investigated the effect of changing viscosity by changing the temperature of 
propylene glycol and glycerine samples. At 50 Hz and 7/p = 35 cm^/sec^, for example, 
we found that in propylene glycol, squares are stable for v < 0.5 cm^/sec while stripes are 
stable for v ^ 0.9 cm^/sec. The crossover region between squares and stripes is complex 
and deserves careful study in its own right. For example, in the same fluid sample, at u 
= 0.63 cm^/sec, we destabilize first to a square state and then, at e = 0.28, to a time- 
dependent, probably chaotic state. At e = 0.34, this state gives way to a stripe state with 
stationary TAM modulations. At e = 0.44, the oscillatory TAM sets in. The wavelength 
of the stationary TAM is 1.4 times the stripe wavelength A, while that of the oscillatory 
TAM is, as mentioned above, 3.3A. Also, by creating a stripe state at high viscosity and 
then slowly decreasing the temperature of the fluid, it is possible to preserve the state 
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even at temperatures where ramping from below threshold leads unambiguously to squares. 
(The reverse is true, as well.) The crossover region is also distinguished by low-frequency 
dynamics for a Od, driven presumably by competition between the two primary states. 
In this regime, we confirmed that stripe states are born with a stationary TAM and then 
further destabilize via oscillatory TAMs at higher e, while square states show 180° phase 
defects. A promising feature of the intermediate-viscosity regime is the possibility, in a 
larger-aspect cell, of finding a spatiotemporally chaotic state at the onset of the primary 
instability [O]. 



IV. DISCUSSION 

Secondary instabilities in viscous fluids that have stripe patterns differ markedly from 
those observed in the low-viscosity, square-pattern case. The TAM instability predicted for 
the latter is in fact only observed for the former. While theory predicts a stationary TAM, 
the instability is oscillatory with a frequency that is generically incommensurate with the 
driving frequency. In the crossover regime, both stationary and oscillatory TAMs are seen. 
It would be interesting to explore the case where ujtam is commensurate with ujr. The closest 
integer ratio, = Sutam, requires higher combinations of u and a than were possible with 
our shaker table. 

An obvious theoretical challenge is to see whether the four coupled amplitude equations 
describing standing waves can account for the oscillatory TAM instability. The coefficients 
in the amplitude equations were derived assuming small damping and small surface slopes 
||13|| . Neither of these conditions is met for high- viscosity fluids although the approx- 
imations are more reasonable for fluids of moderate viscosity. This suggests difficulties for 
perturbation expansions based on a flat interface and small damping, but numerical studies 
remain feasible, as do studies based on models analogous to the Swift-Hohenberg equation 
PJl^ . In the intermediate-damping regime, a small number of amplitude equations may ad- 
equately describe what appear to be spatiotemporally chaotic states. Such a situation has 



recently been discovered in a rotationally anisotropic system |2J], but no isotropic system 
has heretofore been found. 

All of these phenomena — tertiary phase instabilities, locking between primary and 
secondary instabilities, square-stripe competition, the crossover regime — will be much 
easier to explore if one can increase the aspect ratio. Since much larger shaker tables (albeit 
expensive) are available commercially, there is no technical limitation to increasing the aspect 
ratio at least tenfold. 
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FIGURES 



FIG. 1. Patterns observed as the acceleration is ramped slowly from zero. The pictures are 
strobed to show the maximum stripe amplitude, (a) Primary instability at onset. The effect of 
gradients here is significant, (b) Stripe instability, (c) Stripes plus oscillatory TAM instability, 
(d) Phase instability of TAMs. (e) Single defect present in stripe state, (f) Turbulent defects and 
breaking of surface. Glycerine, 44 °C; depth h = 1.0 cm; 7//) = 25 cm^/sec^; u = 1.4 cm^/sec; / 
= 40 Hz, Oc = 6.1 cm^/sec; A = 1.1 cm. F = L/S, where the correlation length ^ is computed using 
the range of linearly unstable wavelengths (see fl^). The estimate of assumes e <C 1, which may 
not be valid for the higher values of e here. 

FIG. 2. Morphology diagram of patterns observed as the acceleration is ramped slowly from 
zero. Regions corresponding to images shown in Fig. |l] are shown. Oil, 22 °C. Similar data was 
obtained using glycerine and propylene glycol. 

FIG. 3. Frequency of the secondary instability. Glycerine, 45 °C, u = 1.2 cm^/sec, 'y/p = 30 
cm^/sec^. The data are also well fit by fx am = a + bf. 
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